CAVENDISH-HEP-2009-23 



Constraining the MSSM with Dark Matter indirect detection data 
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Recently, a claim of possible evidence for Dark Matter in data from the Fermi LAT experiment 
was made by Goodenough and Hooper [8] . We test the Dark Matter properties consistent with their 
claim in terms of the MSSM by a 24-dimensional parameter scan using nested sampling, excluding 
all but a very small region of the MSSM. Although this claim is very preliminary, and not made 
by the Fermi LAT experiment, our scan shows a possible approach for the analysis of future firm 
evidence from an indirect detection experiment, and its potential for heavily constraining models. 



The Fermi Gamma Ray Space Telescope (Fermi), 
launched June 2008 with the Large Area Telescope (LAT) 
experiment on board is one of the most power- 
ful current probes in the indirect search for Dark Mat- 
ter through the detection of its annihilation- or decay- 
products. Gamma ray signatures seem particularly 
generic for Dark Matter models, manifesting as lines from 
directly produced photons, as a hard spectrum from in- 
ternal bremsstrahlung, or a softer spectrum from the sec- 
ondary decays of produced pions. The origin of such 
photons may be cosmological, from the local halo, from 
smaller local structures or from the denser inner region of 
our own galaxy. Fermi data provides an excellent oppor- 
tunity for searches in all these categories [1, 0, 0, 0, 0, 0] • 

Following the public release of Fermi data, Goode- 
nough and Hooper Q have studied the measured gamma 
ray flux from the Galactic Centre, defined to be the re- 
gion 0° < \l\ < 3°, in terms of its angular distribution 
and energy spectrum. They conclude that the data is well 
described by a Dark Matter annihilation scenario, assum- 
ing a Dark Matter particle mass of itidm — 25 — 30 GeV 
and a velocity averaged annihilation cross section of 
(av) ~ 9 x 10~ 26 cm 3 s~\_with the Dark Matter an- 
nihilating dominantly into 66 pairs. 

At this point we would like to express some doubts as 
regards the Dark Matter interpretation of the spectra, 
particularly given that it concerns the notoriously messy 
environs of the Galactic Centre. For the diffuse galactic 
background, Goodenough and Hooper assume what is in 
effect a power-law behavior of the energy spectrum based 
on data from the angular distribution in the region 3° < 
|/| < 6°, and its extrapolation into the Galactic Centre 
as an exponential function. This assumption is clearly 
vulnerable to features in the diffuse background that have 
a stronger dependence on the angular distance to the 
Galactic Centre than the extrapolated function. For a 
recent criticism, see [111 ], The final word on this matter 
will naturally rest with the members of the Fermi LAT 
Collaboration, who are in the best position to understand 
the difficult backgrounds present. 



Despite these doubts, we have found it interesting to 
speculate how the extracted Dark Matter mass and cross 
section can be interpreted in terms of the popular Min- 
imal Supersymmetric Standard Model (MSSM). If firm 
evidence of Dark Matter is indeed found in an indirect 
detection experiment, then features in the spectrum and 
the normalization of excesses point to specific Dark Mat- 
ter masses and annihilation cross sections. It will then 
be natural to ask how this information can be used to 
constrain potential models, such as the MSSM. This is 
the focus of the present Letter. 

Naively, one might expect that a low neutralino mass 
interpretation of the results reported in Q is already ex- 
cluded by collider searches. However, the most commonly 
quoted LEP bound of to^o > 46 GeV @ is based on 
chargino searches and their constraints on the combina- 
tion of low values for fi — M2 , assuming a fixed GUT re- 
lationship between the gaugino masses, Mi w 5M2. The 
invisible width of the Z results in constraints of a similar 
significance on fi and M2, and not directly on the neu- 
tralino mass. In fact, even a nearly massless neutralino is 
not completely excluded in the MSSM, however, a light 
neutralino must be dominantly bino in nature. For a 
recent brief review of the situation, see [Iol |. 

We have used Bayesian sampling techniques to ex- 
plore a 24-dimensional parametrization of the MSSM, as- 
suming diagonal soft supersymmetry breaking sfermion 
mass terms, independent soft mass terms for the gaug- 
inos, non-zero trilinear couplings for the third genera- 
tion sfermions, and keeping /i, tan [3 and mj as free pa- 
rameters in the Higgs sector. We use the nested sam- 
pling Monte Carlo method developed by Skilling [u|, 
with the algorithm implemented in the MultiNest pro- 
gram [3, [lj] . For another use of multi- nested sampling 
in scanning the MSSM parameter space, and a more ex- 
tensive discussion of the method, see [la ]. Our scan re- 
quired approximately one week of continuous running on 
a single modern CPU, though we note that MultiNest 
has a parallel implementation that can be used to gener- 
ate results significantly quicker. 
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The likelihood we use for our scan is constructed from 
the predicted lightest neutralino mass and its velocity 
weighted total cross section, assuming independent mea- 
surements with central values and Gaussian widths of 
m^o = 30±5 GcV and (av) = (9.0±0.9) x 1(T 26 cm 3 s -1 . 
We do not use other observables directly in the likeli- 
hood, besides some direct mass constraints on sparticles 
(see below), as we wish to focus on the effect of the in- 
direct detection observables. For our parameters we use 
flat priors with masses in the range of 10 GeV to 4 TeV, 
with the exception of the gluino mass parameter M 3 for 
which we allow —4 TeV to 4 TeV, the trilinear couplings 
that vary from —1000 to 1000, and tan/3 which lies in 
the interval 2 — 60. 

To evaluate the observables at each given MSSM 
parameter point investigated by MultiNest, we calcu- 
late the sparticle spectrum using the Isasusy patch in 
Isajet 7.78 [Hj], assigning zero likelihood to points de- 
termined to be unphysical. The spectrum is then trans- 
ferred to MicrOMEGAs 2 . 2 [13, EH which is used to calcu- 
late the resulting Dark Matter cross section and density. 

The sampling output consists of posterior PDFs for 
the model parameters, sparticle masses and a few other 
quantities of interest, such as the Dark Matter density. 
In addition, we have implemented the model independent 
LEP direct search bounds on the slepton and chargino 
masses as found in the MicrOMEGAs code, and assigned 
model points in violation a zero likelihood. The resulting 
Xi mass and velocity averaged annihilation cross section 
have Gaussian distributions around the assumed central 
values, thus some a priori allowed models in the MSSM 
are clearly capable of accommodating the input measure- 
ments. Indeed, the point with highest likelihood exactly 
matches the input measurements, and has annihilation 
predominantly to bb pairs. We must now rationalise 
what our most likely models look like, and also determine 
whether they are consistent with other existing particle 
and astrophysical data. 

The marginalised posterior PDF of each weak scale in- 
put parameter is shown in Fig. [T] Of particular interest 
is: i) the tight constraint conferred on M%, due to the 
equally tightly constrained mass and its bino nature, 
and ii) the general preference for higher values of tan/3. 
As expected, the models with the greatest likelihood gen- 
erally show a dominantly bino x° (see the right panel of 
Fig- El), but with a small higgsino component. This is nec- 
essary in order to have a light neutralino that would have 
remained unseen at LEP, which can then only achieve 
the large annihilation cross-section through a Higgs res- 
onance, using 2m^o ~ m^m^t, preferring a small but 
non-zero higgsino component and large tan j3. 

We display the resulting Higgs mass PDFs in Fig. [21 in 
which one indeed observes a resonant annihilation peak 
occurring at roughly twice m^o , split into a double peak 
since the annihilation becomes too efficient when sitting 
exactly on the resonance, leading to a corresponding de- 
crease in likelihood for these models. The Higgs masses 
are also fairly correlated, favouring rrih ttia ~ 2m^o for 



efficient annihilation. The cut off in the rrih distribution 
at high Higgs masses comes from the upper limit on this 
mass, given the squark mass priors, whilst the small rise 
just before the cut off is the result of models in which the 
neutralino acquires a successively greater higgsino com- 
ponent pushing into higher neutralino masses and closer 
to the LEP limits on ii. 

It is evident from Fig. [5] that the LEP bounds on the 
neutral MSSM Higgs masses are amongst the largest hur- 
dles for an MSSM interpretation. Evaluating the precise 
limits given an arbitrary set of weak scale MSSM parame- 
ters is non-trivial, though one can get conservative limits 
by looking at the combined LEP exclusion bounds for 
the various extreme scenarios in 20] . It is interesting to 
note that the weakest bounds at low masses generically 
occur when rrih ~ mA and tan (5 is relatively large, which 
is just what is favoured in our scan. This is due to the 
more difficult to reconstruct Higgs pair-production mech- 
anism, e + e~ — > hA, becoming kinematically allowed and 
competitive with Higgsstrahlung, e + e~ — > hZ. As a re- 
sult we take a very conservative limit of 85 GeV for both 
masses, allowing some points that may be excluded in a 
more rigorous analysis. Even so, this bound clearly ex- 
cludes a large fraction of our models, including the most 
likely. 

Next, we consider the dark matter relic density PDF, 
shown in Fig. [3] (left), which lies significantly lower than 
the value Vih 2 ~ 0.1 indicated by state-of-the-art ACDM 
fits. This is an inevitable consequence of imposing a rea- 
sonably high annihilation cross-section, but does not on 
its own exclude the models. Mechanisms for significant 
boosts in the density through modification of early uni- 
verse cosmology have been suggested, see e.g. [19| , and 
the possibility of multi-component Dark Matter remains. 
Had the models given greater relic densities relative to 
the existing data, one could be more confident in exclud- 
ing them. We also show the PDFs for the neutralino 
components in Fig. [3] (right), confirming its dominant 
bino nature, with some small higgsino admixture and ef- 
fectively no wino component. 

Finally, we consider the impact of the usual selection 
of precision observables, summarised in Table HI Though 
it is difficult to make completely generic statements as to 
why a certain numbers of models fail a certain constraint, 
due to the number of parameters potentially affecting 
each observable, we offer general remarks below. 



TABLE I: Percentage of our posterior samples that fail cuts on 
precision observables. The observables were calculated using 
MicrOMEGAs. 



Observable 


Experimental Limit 


Models excluded 


BR(B -> X s7 ) 
BR(Bs -> nfj) 


(3.55 ±0.42) x 1(T 4 [15] 

< 5.8 x 10" 8 [21] 
(30.2 ±9.2) x 10~ 10 [15] 


46.4% 
41.7% 
99.7% 



The branching ratio for the decay B — > X s j is 
tightly constrained experimentally to lie close to the SM 




FIG. 1: Posterior PDFs for all weak scale input parameters. 
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FIG. 2: Posterior PDFs for m h (left) and m A (right). 




FIG. 3: Posterior PDF for Dark Matter density Q.h 2 (left) 
and neutralino make-up (right). 



value, but can be modified within the allowed range by 
SUSY contributions, predominantly involving loops with 
charged Higgs bosons and charginos. The chargino con- 
tribution is enhanced at large tan [3 and has potential 
large logarithms of msusy Jrnw- For the allowed range 
we use the evaluation in [15j |. A significant number of 
our models give BR(B — > X s ^) values higher than the 
measured upper limit, in keeping with the observed pref- 
erence for higher tan/?. 

Similarly, the branching ratio for the decay B s — » fifi 
can be enhanced in the MSSM via interactions involving 
neutral Higgs bosons. The effects are again higher at 
large tan /?, and one would expect an increase with lower 
neutral Higgs masses. A large number of our models 
are excluded by the the CDF Run II constraint on this 
branching ratio [2l|, including our best fit model, and 
those that remain generally have lower values of tan 0. 

By far the largest impact on our models comes from re- 
quiring the correction to the muon anomalous magnetic 
moment, 5a M , to account for the observed discrepancy 
with the SM. Almost all of our posterior points are in- 
consistent due to having too low Sa^ values. The main 
SUSY contributions to a M result from smuon-neutralino 
and sneutrino-chargino loops (22J. Our models generally 
prefer heavy neutralinos and charginos Xii f° r i > U due 
to a large M2, pushed to higher values to maintain the 
specific character of the Xi required by our input data, 
while the limit on the lightest chargino mass prevents fi 



from becoming very low. As a check we have taken our 
highest likelihood model, lowering the smuon mass pa- 
rameters whilst also lowering M2 , which indeed increases 
Sa M . 

If one discounts the apparent 3cr discrepancy between 
the experimentally observed a M and the SM prediction, 
the posterior survives more or less unscathed. A thor- 
ough investigation of each of these precision observables 
would no doubt reveal further sensitivities and correla- 
tions but we considered this beyond the scope of this 
short Letter. 

Given the above, it is reasonable to ask if any MSSM 
models consistent with all current experimental limits 
on supersymmetry remain compatible with the Goode- 
nough and Hooper claim. Indeed 31 points, out of a to- 
tal of 45301 posterior points, survive after imposing the 
bounds in Table [I] and in fact all have reasonable likeli- 
hoods, having annihilation cross-sections and neutralino 
masses close to the input values. These models have a 
significantly lower mass for the first and second genera- 
tion sleptons, as one would expect from imposing the Sa^ 
constraint. The most likely point remaining after impos- 
ing the conservative LEP Higgs bound detailed above 
has ttidm — 31.8 GeV, a velocity averaged annihilation 
cross section of (av) ~ 8.73 x 10~ 26 cm 3 s _1 and Higgs 
masses (m^myi) = (96.8,96.5) GeV, and is thus not too 
far from satisfying the Goodenough and Hooper claim. 
There are a few models with Higgs masses greater than 
110 GeV, but these are considerably less likely. If one 
ignores the lower bound on ^a^, the most likely model 
has m DM 01 30.3 GeV and (av) ~ 9.04 x 10" 26 , and 
(m h ,m A ) = (106.9,106.5) GeV. 

In conclusion, we have shown that the Dark Matter 
interpretation of the Fermi LAT data from the Galac- 
tic Centre in Q is generally hard to reconcile with neu- 
tralino Dark Matter in the MSSM. However, a small re- 
gion of parameter space remains with a very light bino 
neutralino with some higgsino admixture, Higgs masses 
mh and near the border of the LEP excluded re- 
gion, fairly light sleptons and moderate values of tan/3. 
This demonstrates the power of multi-modal nested sam- 
pling — using only approximate information on the Dark 
Matter mass and annihilation cross section, as expected 
from a indirect detection signal — to constrain mod- 
els with a large number of parameters and complicated 
degeneracies, performing an efficient scan over a many- 
dimensional parameter space. 

Our current scan was very efficient in restricting the 
MSSM parameter space due to the difficulty in the MSSM 
of simultaneously getting a low neutralino mass and a 
large annihilation cross section, while at the same time 
evading earlier collider bounds. We therefore feel that 
one would do well to repeat this interesting case study if 
more definite evidence of Dark Matter is found in indirect 
detection experiments, taking care to also incorporate 
properly the astrophysical uncertainties associated with 
the halo distribution. Finally, although it is possible to 
get close to accommodating the Goodenough and Hooper 
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claim for a small fraction of our models, dedicated fans 
of neutralino Dark Matter in the MSSM should probably 
hope that the present interpretation is a false alarm. 
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